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We present the first measurement of mid-rapidity <p vector meson production in 
Au Au collisions at RHIC GeV) from the STAR detector. For 

the 11% most central collisions, the slope parameter from an exponential fit to 
the transverse mass distribution is T= 379±51(stat)± 45(syst) MeV, the yield 
dN/dy = 5.73±0.37(stat)± 0.57(syst) per event and the ratio is found to 

be 0.021±0.001(stat)± 0.004(syst). We currently place the value of the N,f,/Ni^- 
ratio between 0.10 and 0.16. The measured ratios N^/Nf,- and N^/Nj^- , as well 
as T for the 0 meson at mid-rapidity do not change for the selected centrality bins. 


1 Introduction 


The central topic of relativistic heavy ion physics is the study of Quantum 
Chromodynamics (QCD) in extreroe conditions of high temperature and high 
energy density over large volumes Jj. Vector mesons may probe the dynamics 
of particles and chiral symmetry □ in relativistic heavy ion collisions: their 
production mechanisms and suhseouent dynamical evolution have been a topic 
of experimental investigation HuRQ. The </> meson is of particular interest due 
to its ss valence quark content, which may make sensitive to strangeness 
production from a possible early partonic phase Lrffl. 

In central Pb-\-Pb collisions at the CERN SPS (nucleon-nucleon center of 
mass energy ~ 17 GeV), the slope parameter (T) in an exponential fit 

to the transverse mass (mt) distribution at mid-rapidity (oc follow^ 

systematic trend as a function of hadron mass for pions, kaons and protons td. 
This observation is indicative of a common expansion velocity developed in 
the final state for pions, kaons, and protons £ 3 . Thej-slope parametep, how¬ 
ever, measured for multi-strange hyperons S and Ul £ 3 , and for J/V’lij show 
deviations from a linear mass dependence, suggesting that these particles do 
not interact as strongly in the final state at SPS en^gies £ 3 . Measurements of 
(p meson production at the SPS were inconclusive u’cl. Significantly different 
values for the (p slope parameter have been obtained from exponential fits to 
the measured rat spectra in central Pb-\- Pb collisions wheruusing the K^K~ 
decay channel a and when using the ll~ decay channel □ of the_di_peson. 
This difference, however, is not apparent in peripheral collisions £ 3 £ 3 . Pos- 

* scenarios to explain the difference have been discussed in the literature 




We report the first measurement of mid-rapidity (|y| < 0.5) (j) production 
in Au + Au collisions at RHIC i^/s^ = 130 GeV) via the (j) K~ decay 
channel (branchir^ ratio = 0.491) using the Solenoidal Tracker At RHIC 
(STAR) detector E3. Systematics of (j) meson production as a function of 
centrality at RHIC as well as its dependence will be discussed. 

2 Analysis 

The STAR detector consists of several detector sub-systems in a large 
solenoidal analyzing magnet. For the data taken in the year 2000 and pre¬ 
sented here, the experimental setup consisted of a Time Projection Chamber 
(TPC), a Central Trigger Barrel (CTB), and two Zero Degree Calorimeters 
(ZDC) located upstream along the beam axis. The TPC is a cylindrical drift 
chamber with multi-wire proportional chamber readout. With its axis aligned 
along the beam direction, the TPC provided complete azimuthal coverage. 
Surrounding the TPC was the CTB, which measured energy deposition from 
charged particles. The ZDC’s measured beam-like neutrons from the frag¬ 
mentation of colliding nuclei. The CTB was used in conjunction with the 
ZDC’s as the experimental trigger. 

Data used in this analysis were taken with two different trigger conditions: 
a minimum-bias trigger requiring a coincidence between both ZDC’s and a 
central trigger additionally requiring a high hit multiplicity in the CTB. The 
central trigger corresponded to approximately the top 15% of the measured 
cross section for Au + Au collisions. Data from both the minimum-bias trigger 
and central trigger were used for this analysis. 

Reconstruction of the (f> was accomplished by calculating the invariant 
mass (rriinv), transverse momentum (pt), and rapidity (y) of all permutations 
of candidate K^K~ pairs. The resulting rriinv distribution consisted of the 4> 
signal as well as combinatorial background. The shape of the-jCombinatorial 
background was calculated using the mixed-event technique 

For the centrality measurement, the raw total charged multiplicity distri¬ 
bution within a pseudo-rapidity window |? 7 | < 0.75 was divided into three bins 
corresponding to 85-26%, 26^^ and the top 11% of the measured cross sec¬ 
tion for Au + Au collisions Events were selected with a primary vertex 

z position (z) from the center of the TPC of \z\ < 80 cm. These events were 
further divided according to z in 16 bins, and event mixing was performed for 
events within each bin to construct background distributions with reduced 
acceptance-induced distortions in the mixed-event background. Consistent 
results were obtained when we constructed the background distribution using 
like-sign pairs from the same event. 
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Figure 1. Left: The measured {dEjdx —> vs. p for reconstructed tracks in the TPC. 
Right: Reconstruction efficiency vs. mt — rriff, vs. rapidity for minimum-bias events. 


Particle identification (PID) was achieved by correlating the ionization 
energy loss {dE/dx) of charged particles in the TPC gas with their measured 
momentum (left panel of Figure ||). By truncating the largest 30% dE/dx 
values along the track, a sample was selected to calculate the mean {dE/dx). 
For the most central events, the average (dE/dx) resolution was found to be 
about 11%. The measured (dE/dx) is reasonably described by the Bethe- 
Bloch function smeared with a resolution of width a. Tracks within 2a of the 
kaon Bethe-Bloch curve were selected for this analysis. 

To obtain the (j) spectra, same event and mixed event distributions were 
accumulated and background subtraction was done in each (mt,y) bin. The 
mixed event background rriinv distribution was normalized to the same event 
rriinv distribution in the region above the (/> mass (1.04 < rriinv < 1-2 GeV/c^). 
The raw yield in each bin was then determined by fitting the background sub¬ 
tracted rriinv distribution to a Breit-Wigner function plus a linear background 
in a limited mass range (left panel of Figure A slight mismatch between 
the same-event and mixed-event background distributions, due to Coulomb 
interactions, track merging and otheuesidual correlations leads to a structure 
in the subtracted mass distribution EH. The width of the fit to the invariant 
mass distribution is consistent with the natural width of the (j) convoluted 
with the resolution of the TPC. 

The resulting raw (j) yield for each mt, y and centrality bin was then cor¬ 
rected for tracking efficiency and acceptance using Monte Carlo simulations 
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Figure 2. Left: Invariant mass distribution for candidate K'^K~ pairs after background 
subtraction for the 11% most central collisions. The width of the invariant mass mass 
distribution is consistent with the natural width of the cj) convoluted with the resolution of 
the TPC. Right: Transverse mass distributions of cj) from Au + Au collision at = 130 

GeV for three centrality bins. Dashed lines are exponential fits to the data. For clarity, data 
points from the 11-26% and 26—85% centrality bins are scaled by 0.1 and 0.01, respectively. 
Error bars shown are statistical only. 


of physics processes and detector response. The reconstruction efficiency de- 
pended on the pt of the (j)^ ranging from ~ 10% at pt = 0.46 GeV/c up to 
^ 40% at Pt = 1.4 GeV/c (right panel of Figure The PID efficiency cor¬ 
rection for the (p was calculated as the square of the single kaon PID efficiency 
and included the centrality dependence of the dE/dx resolution. The cor¬ 
rected (j) invariant yields for three centrality bins are shown in the right panel 
of Figure All results presented here are for reconstructed (p mesons within 
one unit of rapidity centered around y = 0 (|y| < 0.5) and 0.46 < pt < li4 
GeV/c. In the region where the pion band crosses the kaon band in dE/dx c3, 
corresponding to the kaon pt ~ 0.8 GeV/c, the signal to background ratio de¬ 
grades. This leads to the larger statistical error bars in the most central bin 
and prevented the extraction of the p yields in this region for the two lower 
centrality bins. The spectra were fit to an exponential 


1 d'^N 
27rmt drritdy 


dN/dy 

27rT(m0 -I- T) 


( 1 ) 


with the slope parameter T and yield dN/dy set as free parameters. The 
results obtained are listed in Table |l|. The fraction of p mesons in the measured 







Event Centrality 

0-11% 

11 26% 

26-85% 

T(MeV) 


379±51 

369±73 

417±76 

dN/dy 


5.73T0.37 

3.33T0.38 

0.98T0.12 

N^/Nh- 


0.021T0.001 

0.019T0.002 

0.019T0.002 

N^/Nk- 


4= 

0.10 TO 0.16 


All numbers are preliminary 


Table 1. Mid-rapidity 0 slope parameters T, extrapolated yield dN/dy and the ratios 
and for three centrality bins. Both the statistical and systematic un¬ 

certainties are listed. 


Pt region assuming an exponential distribution is ^ 70%. Also l^ed is the 
midrapidity ratio of the </> yield to the negative hadron {h~) yield E3 for three 
multiplicity bins. 

The major systematic uncertainties for this analysis include contributions 
from PID efficiency and tracking efficiency. By varying PID and track quality 
requirements, we estimate a systematic uncertainty of ± 12% for T and ± 
10% for dN/dy. Systematic errors for the N^/Nf^- ratio also includes the 
uncertainty in the h~ yield. The full range of the systematic uncertainty for 
N^/Nh- is ±20%. 


3 Results 

For the most c^tral heavy ion collisions, there is an increase in T from 
the AGS H to SPS 0 to RHIC. The slope parameters, however, from p + p 
colliskps show no significant dependence on collision energy up to ^/s = 63 

GeV aid. 

Since the (j) and ^ti-proton have similar masses and very different scat¬ 
tering cross sections B, comparison of the spectraf shapes would shed light 
on collision dynamics. In the most central Au + Au collisions at RHIG, the 
4> slope parameter is 379±51(stat)± 45(syst) and there is no dependence on 
event centrality (Table within our statistical uncertainty. The anti-proton 
slope parameter using the same fit function, however, measured in the pt 
range 0.25 < pt < 1 GeV/c and without correction for feed-down from anti- 
hyperons, is found to be ovei 150 MeV higher than the (j) meson slope measured 
in 0.5 < Pt < 1.7 GeV/c M. Note that if a strong collective flow develops in 
the system, the measured slope parameter should depend strongly on the fit¬ 
ting range. Measurements of the (j) and anti-proton over a much broader range 
of Pt will yield a more definitive conclusion on the dynamics of these particles. 
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Figure 3. \/Snn dependence of the mid-rapidity Left: N(j)/Nj^- and Right: Ntp/Nf^- 
ratios. The data points are from the ^ K~ decay channel. Filled symbols represent 

the results extracted from the most central heavy ion collisions and the open symbols 
represent the results from p-\-p (17 and 63 GeV points) and p-\-p (1800 GeV) collisions. 
The shaded region in the plot of the ratio represents the estimate of the ratio in 

= 130 GeV Au-\-Au collisions at RHIC. Error bars shown are statistical errors only. 


The energy dependence of the mid-rapidity ratio is shown in the 

left panel of Figure |[ The shaded region corresponds to our best estimate 
of the N^/Nk- ratio in = 130 GeV Au + Au collisions at RHIC. The 

N^/Nr- ratio has no measurable dependence on collision energy and system 
size over three orders of magnitude from Au + Au collisions at the AGS Q to 
p + p collisions at the Tevatron Ed. 

The energy dependence of the mid-rapidity N(j,/Nf^- ratio is shown in the 
right panel of Figure ^ In heavy ion collisions, increases with colli¬ 

sion energy indicating that (j) production increases faster than h~ production 
up to = 130 GeV. Although there seems to be a sigm^ant increase in 

N^/Nfi- ratio from p-bp collisions between 17 and 63 GeVB’td, the statistical 
uncertainty in the 63 GeV point is too large to determine the energy depen¬ 
dence. Note that the ratio ^ Tevatron energies (p -b p at a/s = 1800 GeV) 
was found to be about 0.01 Ed. 










4 Conclusion 


In summary, using the STAR detector we have measured mid-rapidity </> 
production from Au + Au collisions at ^/s^ = 130 GeV. In the most central 
collisions, the (j) slope parameter, T= 379±51(stat)± 45(syst) MeV, is lower 
than that of anti-protons in the measured pt region. Within statistical un¬ 
certainty, there is no variation in (j) slope parameters and the ratio N^/Nf^~ 
for the selected centrality bins. There is no measurable dependence of the 
Nip/Nr- on collision energy and system size spanning three orders of magni¬ 
tude. The (j) slope parameter and the ratio N^/N^- increase from ~ 5 

to 130 GeV. 
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